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ABSTRACT 

Aims. Solar radio type II bursts are rarely seen at frequencies higher than a few hundred MHz. Since metric type II bursts are thought 
to be signatures of propagating shock waves, it is of interest to know how these shocks, and the type II bursts, are formed. In particular, 
how are high-frequency, fragmented type II bursts created? Are there differences in shock acceleration or in the surrounding medium 
that could explain the differences to the "typical" metric type lis? 

Methods. We analyse one unusual metric type II event in detail, with comparison to white-light, EUV, and X-ray observations. As the 
radio event was associated with a flare and a coronal mass ejection (CME), we investigate their connection. We then utilize numerical 
MHD simulations to study the shock structure induced by an erupting CME in a model corona including dense loops. 
Results. Our simulations show that the fragmented part of the type II burst can be formed when a coronal shock driven by a mass 
ejection passes through a system of dense loops overlying the active region. To produce fragmented emission, the conditions for 
plasma emission have to be more favourable inside the loop than in the interloop area. The obvious hypothesis, consistent with our 
simulation model, is that the shock strength decreases significantly in the space between the denser loops. The later, more typical type 
II burst appears when the shock exits the dense loop system and finally, outside the active region, the type II burst dies out when the 
changing geometry no longer favours the electron shock-acceleration. 
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1. Introduction 

Radio type II bursts are observed in association with flares and 
coronal mass ejections (CMEs). These bursts can be observed 
at metric wavelengths (coronal type II bursts) and at decame- 
ter and longer wavelengths (interplanetary type II bursts). The 
mechanism behind the bursts is generally assumed to be a propa- 
gating shock which creates electron beams that excite Langmuir 
waves, which in turn convert into radio waves at the local plasma 
fre quency and its second harmonic ( [Melrose 1980t [Cairns et 



al. 20Q3| ). As shocks can be formed in various ways (for an 



overview see, e.g., Warmuth, 2007; for the terminology see, e.g., 
Vrsnak, 12005), it is not evident that all solar radio type II bursts 
are formed in the same way. 

Most of the interplanetary type II bursts are thought to be 
created by CME-driven shocks, but there is observational ev- 
idence that at least some coronal metric type II bursts are ig- 
nited by smaller-scale processes associated with the flare en- 
ergy release, such as high-speed plasma jets/blobs or loop ejec- 
tions/expansions ("Klein et aL| '1999'; "Pohjolainen et al. 



Khan & Aurass , 2002(]Klassen et al., 2003 ; Daup hin et al^ 



2001 



2006 



Pohjolainen 2008| ). Statistically, there is support for the idea that 



metric ty pe II bursts have the ir root cause in fast coronal mass 
ejections ([ Cliver et"aL| | 1999 |), but also th at they are not caused 
by shocks driven in front of CMEs ( Cane & Erickson[|2005| . 

Metric type II bursts can be observed in dynamic rad io spec- 
tra as slowly (0.1 - 1.0 MHz s"^) drifting emission lanes ( [Nelson 
[& Melros"e|[1985 )). The start frequency of metric type II bursts 
is usually at about 100-200 MHz, and the bursts have a typi- 



cal duration of 2- 18 min ( [Subramanian & Ebenezer|[2006[ ). In 
some cases it is possible to separate a slowly drifting "backbone" 
in the burst emission, with fast-drifting (-10 MHz s"^) emis- 
sion stripes shooting up and down from the backbone ( Mann & 
Klassen, 2005). These features have been named as "herring- 
bones", and they bear some resemblance to radio type III bursts 
which are signatures of outflowing electron beams. However, the 
drift rates of type III bursts hav e been found t o be higher than 
those of herringbones ( Mann & Klassenl|2002|). 

In this paper, we analyse in detail one metric type II burst 
that occurred on 13 May 2001. The burst started at an unusu- 
ally high frequency, proceeded showing fragmented and curved 
emission bands, which were visible at the fundamental and sec- 
ond harmonic plasma frequencies. As the radio burst was associ- 
ated with a flare and a CME, we do a multi- wavelength analysis 
in order to identify the origin of the shock that formed the radio 
type II burst. We then use numerical MHD simulations to verify 
the results from the data analysis. 



2. Observations 

A GOES M3.6 class flare was observed to start at 02:58 UT on 
13 May 2001 in NOAA AR 9455 (S18W01). The flare max- 
imum occurred at 03:04 UT, and at 03:50 UT a coronal mass 
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ejection was detected in the LASCO (Brueckner et al. 1995 ) C2 
images. The CME front moved toward the South at a speed of 
430 km s"^ (LASCO CME Catalog, both the linear and second- 
order fits to the height-time data give similar speeds). Backward- 
extrapolation of the heights indicates a CME start time near 
03:00 UT. Fig.[T]shows the GOES soft X-ray flux curve, with the 
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Fig. 1. Left: GOES soft X-ray flux curve at 02:50-04:20 UT, with flare class classification A-X. The dashed region indicates 
Yohkoh satellite night time, when no observations are available from the experiment. Plane-of-the-sky heliocentric CME front 
distances (boxes combined with a dotted line) are from the LASCO CME Catalog. Black boxes between 03:01:30 and 03:03:30 
UT show the projected distances of the soft X-ray loop-like front, measured from the center of the eruption region at 1 R©. Right: 
Detailed plot of the estimated distances/heights measured in the time interval 03:01:00-03:04:20 UT. Projected distances (km) 
from the eruption center are shown for the SXT loop-like front (stars/dash-dotted line) and the EUV 'blob' (crosses/dotted line). 
The estimated type II burst heights are calculated with the lOxSaito density model for the fragmented part (the three height- 
time data points, marked with diamonds, are from the start times of each three fragment. The type II heights for the later regular 
part (dashed lines) are calculated with the hybrid density model (corresponding to SxSaito) which describes a "standard" solar 
atmosphere. Applying a high-density model like lOxSaito for the regular type II part at < 130 MHz would imply at least streamer 
densities, see text for details. Yohkoh HXT hard X-ray counts in the 23-33 keV energy channel are also shown (bottom panel). 
Note that the drop in counts near 03:04 UT is due to the satellite entering night time. 



CME plane-of-the-sky heigh ts during 02:50-04:20 UT. Yohkoh 
HXT ( [Kosugi et al.[ |1991| ) hard X-ray counts are shown at 
03:01:00-03:04:20 UT, before satellite night ended the obser- 
vations. 

Radio emission at decimetric-metric wavelengths was ob- 
served to start at 03:01:45 UT. It consisted of a type II burst 
that was clearly visible as two emission lanes, starting near 500 
and 1000 MHz. It was followed by bursty continuum emission 
at 700-300 MHz. These features are shown in the HiRAS dy- 
namic spectrum that covers a wide, 2500-25 MHz frequency 
range, in Fig. [2] The type II burst emission ended abruptly at 
03:05:15 UT near 90 MHz at the fundamental emission band 
(visible also in the RSTN observations at 180-25 MHz). No in- 
terplanetary type II emission was o bserved at lower freq uencies, 
as verified from the Wind WAVES ( [Bougeret et al. 1995| ) obser- 
vations at 14 MHz-20 kHz. 



2.1. Fragmented radio type II burst 

The radio type II burst showed a pair of fundamental and second 
harmonic emissions, which clearly marks it as plasma emission. 



HiRAS 



fp = 9000 ^/J^e, 



(1) 



where the local electron plasma frequency fp is expressed in Hz 
and the electron density rie in cm"^. As plasma emission at the 
fundamental is directly related to the local electron density, we 
can use the spectral observations to calculate electron densities 
and to estimate burst source heights and speeds. Height esti- 
mates require, however, the use of atmospheric density models, 
which are not unambigu ous, see e.g., |Vrsnak et al.| ( [2004aj and 
Pohjolainen eFS?! pOD?] ). 

The 500 MHz start frequency (fundame ntal plasm a emis- 
sion) is unusually high for a type II burst ( Lin et alT] [2006). 
During the first two minutes of the burst, the emission was frag- 
mented into separate, curved emission bands. The frequency 
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Fig. 2. HiRAS dynamic spectrum from the 25-2500 MHz fre- 
quency range, at 03:00-03:10 UT on 13 May 2001. 'F' notes 
emission at the fundamental and 'H' at the second harmonic 
plasma frequency. Arrows point to the fragmented emission 
bands and dashed white lines outline the later- appearing "reg- 
ular" type II burst lanes. 



drifts (fundamental emission) within the fragments were be- 
tween 1.8 and 4.3 MHz s"^ These drift rates are high com- 
pared to type II bursts in general, but consistent with drift rates 
found for high-frequency type II bursts, see Fig. 2a in Vrsnak] 



|et al. p002 ). The first "fragment" showed emission between 



500 and 420 MHz, which correspond to densities in the range 
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of 2-3 X 10^ cm ^. The second fragment near 03:02 UT gives 
densities 1 - 2 x 10^ cm"^ (400-310 MHz), and the third near 
03:03 UT gives densities 2-6 x 10^ cm"^ (220-130 MHz). 
These values indicate that the source regions were dense, similar 
to active region loops. Rough estimates for heliocentric heights, 
using atmospheric density models up to ten-times Saito ( Saito[ 
1970), give values in the range of 1.03-1.43 R©. These esti- 
mates are suggestive, since the models do not give heights for 
erupting structures. The calculated type II burst source heights 
(lOxSaito densities) are shown at the start times of the three 
fragmented emission bands in Fig. [T] The burst velocity, using 
these heights, is as high as 2 300 km s"\ which probably re- 
flects the gradual density gradient in the Saito density model. 
For comparison, if we calculate the burst speed from a steeper 
density gradient, using the emission frequency (500 MHz) and 
the general drift-rate of the bands (2.4 MHz s"^), and the scale 
height of 73 000 km (hydrostatic isothermal scale height at 1.25 
MK coronal temperature and heliocentric height of 1.1 R©), we 
get a burst speed of 1 000 km s"^ (for the method see, e.g., 
Pohjolainen et al., |2007|| ). 

At 03:03:35 UT the type II lanes got wider and started to 
look like a typical type II burst. The start frequency of this "reg- 
ular" part was 130 MHz at the fundamental, and the frequency 
drift was about 0.4 MHz s"^ This is within the usual drift rates 
observed for type II bursts (between 0.1 and 1.0 MHz s"^). To 
estimate the height for this part of the type I I burst we use the 
hybrid density model by |Vrsnak et al.| ( |2504a ). This model gives 
density values similar to five-times Saito in the low corona but 
also links coronal densities to those in the interplanetary space. 
The heliocentric heights for the burst source are 1.22 R© at the 
beginning of the burst (130 MHz), and 1.34 R© at the end of 
the burst at 03:05: 15 UT (90 MHz). A height-time trajectory for 
this regular type II burst part is also shown in Fig. [T] Applying 
a high-density model like lOxSaito also for the regular type II 
part would make the two type II burst tracks converge, but this 
would imply that the regular part was formed at least at streamer 
densities, above 1.4 R©. If the burst source was moving along the 
density gradient the deduced speed from the hybrid model is ap- 
proximately 840 km s"^ If the direction of motion diff'ered from 
that, the speed could have been higher (an angle of 45° between 
the driver velocity vector and the density gradient produces a 
speed 1.4 times higher). 



2.2. Eruptive flare 



100 



The pre-eruption Yohkoh SXT (Tsuneta et al. 1991 ) images 



show several loop systems, connected both to the active region 
and regions nearby (Fig. [3|. Small-scale pre-flare brightenings 
were observed before the flare start at 02:58 UT. The flare re- 
gion is indicated with an arrow in Fig.|3]_ 



EUV images taken with TRACE ( Handy et al.| [19991 ) show 
a filament eruption, in which most of the material is moving to- 
ward the Southeast. The outermost front of the filament moves 
with a projected speed of about 380 km s"\ but there is also a 
separate 'blob' that moves more to the Southwest with a pro- 
jected speed of 450 km s"\ see Fig. |4] The TRACE image in 
Fig.|4]is at 03:02:04 UT, very near the time when the first frag- 
mented radio type II burst band appeared. 

Yohkoh SXT observed the flare until 03:04 UT, when satel- 
lite night set in. The partial SXT frames between 03:02 and 
03:03 UT reveal a loop-like front moving Southward at a pro- 
jected speed of about 650 km s~, see Figs. and m The pro- 
jected distances for the SXT loop-like front and the EUV blob, 
calculated from the eruption center where loop footpoints were 
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Fig. 3. Yohkoh SXT composite of images taken with two difl'er- 
ent filters and field of views at 02:54 UT on 13 May 2001. It 
shows the pre-eruption flare region and interconnecting loops. 
Axes are in arc seconds and white areas in the image mark satu- 
rated pixels. 
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Fig. 4. TRACE EUV image at 03:02:04 UT (top), showing the 
erupting filament. Most of the filament material moved toward 
the Southeast (speed 380 km s"^), but a blob-like structure 
moved more to the Southwest (speed ^ 450 km s"^). Arrows give 
the direction of motion for these structures. Yohkoh SXT image 
(bottom) shows the same region in soft X-rays. The pointing for 
TRACE seems to be off' from that of Yohkoh and so point-to 
point comparison is ambiguous. 



visible prior to the event, are shown in Fig. [T] together with 
the type II burst height estimates (at lOxSaito atmospheric den- 
sities). The height of the first type II burst fragment is low 
compared to the SXT and EUV structure heights, but there is 
also some uncertainty in the true density gradient. The frag- 
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Fig. 5. SXT difference images (AlMg filters) at 03:02:50 UT and 
03:03:30 UT, showing the loop-like eruption front in soft X-rays. 
The later difference image is near the time of the appearance of 
the 'regular' metric type II burst at 03:03:35 UT. The projected 
plane-of-the sky speed of the soft X-ray front is about 650 km 
s"^ (arrows point to the loop front). 

ment heights, if true, suggest that the shock may have travelled 
through the EUV and/or soft X-ray structures. On the other hand, 
a low type II start height also suggests that the type II burst may 
have been emitted at the low-lying flanks of the shock. If, how- 
ever, the speed of the propagating shock remained in the range of 
1 000 - 850 km s"^ during the whole event, the fragmented burst 
could have started at heliocentric height of 1 . 1 R© (^^ 70 000 km 
above the photosphere), which is just above the SXT loop-like 
front. 

As Yohkoh SXT was observing with two different filters. 
Be and AlMg, it is possible to es timate the emission me asure 
EM using the filter ratio, see, e.g.,|McTiernan et aL](|1993|). The 



emission measure for the outermost SXT loop-like structure at 
03:02:18 UT is 6 x lO^^ - 10^4 cm"^ per pixel. As one pixel is 
2.5'' and the loop width is approximately 10'' in the images, we 
can estimate the loop depth along the line of sight to be equal 
to the loop width, ^ 4 pixels, and thus get an emitting volume 
y of 2.3 X 10^^ cm^ per pixel. Assuming an isothermal homo- 
geneous source, the electron density N can be calculated from 
N ^ ^|EM|V, which gives a value of about 1.6-2.1 x 10^ cm"^ . 
For comparison, the fragmented radio emission near 03:02 UT 
indicated electron densities of 1 - 2 x 10^ cm"^, a value similar 
to our soft X-ray loop densities. 

The SXT diff'erence image at 03:03:30-03:02:50 UT 
(Fig.[5 ) shows that the (projected) soft X-ray loop-like front was 
located approximately 1 80" from the center of the flaring region. 
This corresponds roughly to a distance of 126 000 km. For com- 
parison, the estimated height of the later, regular part of the type 
II burst source at 03:03:35 UT is about 152 000 km. Even taking 
into account projection effects, the radio burst source is located 
higher than the soft X-ray loop-like front at that time. 



3. Coronal conditions 

If we assume a structure moving through a system of dense loops 
at a speed of 850 km s"^ (similar to the inferred metric type II 
shock speed; see also Sect. 5), shock- acceleration could be cre- 
ated if the speed of the disturbance exceeded the local magne- 
tosonic speed. In the solar corona, the local Alfven speed va is a 
good approximation for the magnetosonic speed. 



2x lO^i^ 



(2) 



locity super- Alfvenic at the density of the first type II fragment 
(2-3x10^ cm"^), the magnetic field strength should be less than 
22 G. At the time of the third fragment, the electron density was 
2-6 X 10^ cm"^, implying a field strength less than 10 G. The 
empirical scaling law presented by D ulk & McLean, ( jl978j ) gives 
the strength of the coronal magnetic field above active regions as 
B - 0.5 xh~^-^ G, where h is the height above the photosphere 
expressed in solar radiae. Field strengths of 22 G would then 
appear at about heights of 56 000 km and 10 G at about 95 000 
km. The outermost SXT loop heights, approximated from the 
observed projected distances between the loop footpoints and 
the loop fronts, were about 66000 km at the time of the first 
type II fragment and about 100000 km at the time of the third 
fragment. The EUV blob was located lower, at about 47 000 and 
7 1 000 km, respectively. The estimated field strengths imply that 
the shock could have been located at about the SXT loop heights 
at the time of the fragmented type II burst. 

Basically, any propagating shock would be able to excite 
a type II burst, since the type II characteristics mainly reflect 
the medium where the shock is propagating. The SXT loop-like 
front could have been a low-coronal signature of the CME fRust 
|& Hildnerl \T916\ [Vrsnak et al.| |2004b| [Temmer et al.| |2Q08| ), 
which drives the shock. The shock could also have been ignited 
by the flare, and propagated as a blast wave. In this scenario, if 
the SXT loop-like front roughly corresponds to the CME front, 
the type II fragmented emission would be excited at segments 
where the shock overtakes some part of the CME - and the type 
II fragments would be attributed to the inhomogeneities within 
the erupting structure. One difference between shock types is 
that a blast wave should weaken in time as it does not gain energy 
on the way, and because it also expands the energy is divided 
into a larger area. Eventually a blast wave shock will become so 
weak that it cannot accelerate electrons any more, leading to a 
slow fade out of the type II burst. 

The type II burst disappearance from the spectrum was 
quite sudden at 03:05:15 UT. The abrupt end suggests that the 
shock may not have been weakening but instead it entered a 
region where the magnetic field orientation was very differ- 
ent, e.g., there was a transit from a perpendicular to a parallel 
shock region. Shock-drift acceleration is very efficient for quasi- 
perpendicular shocks, and a change in the geometry could stop 
accelerating particles. Of course, this could happen both to a 
driven shock and a freely propagating blast wave shock (before 
it weakens). In the case where a blast wave shock overtakes a 



CME, type II emission could end after the overtake ( [Gary et aL 
.1984J . 

4. MHD modeling 

We utilize numerical simulations to study the shock structure 
induced by an erupting CME in a model corona including dense 
loops. Since we are interested in modeling the first minutes of 
the eruption, we consider a local model of the low corona. The 
time-dependent, ideal MHD equations augmented with gravity 
are solved in a two-dimensional Ca rtesian domain. The d etails 
of the model have been presented bylPomoell et al. (2008 ), who 
employed a largely similar model to study the waves excited by 
an erupting flux rope. 

The magnetic field configuration consists of a flux rope and 
an arcade-like background field. The field of the detached flux 
rope is created by a line current of radius ro situated at distance 
h above the solar surface, and the arcade structure of the back- 



where B is the magnetic field strength and Ug is the local electron 
density, in gaussian units. In order to make the disturbance ve- 



gro und st ructure follows the coronal arcade model of [Oliver et 
al. (1993]). While the density of the ambient corona decreases 
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Density Speed 



Table 1. Alfven speed in km/s at four locations (x, y) for the two 
initial conditions 



Run 


(-0.16,1) 


(0,1.1) 


(0,1.2) 


(0,1.45) 


A 


335 


515 


220 


571 


B&C 


457 


703 


30 


779 



exponentially due to the choice of constant gravitational accel- 
eration, gQ = 2.74 X 10"^ cm s"^, and constant initial temperature 
To = 1.5 X 10^ K, the density inside the flux rope is enhanced to 
model denser filament material. 

To include coronal loops with an increased density into the 
model, we apply a simple approach in which the magnetic field 
strength is decreased in the denser loops so that the total pressure 
is constant at a constant height . For the details of th is procedure, 
we refer to the work of Odstrcil & Karlicky|pOOO|). Note that in 
real coronal loops, both the magnetic field strength and density 
are enhanced compared with the ambient corona. We discuss in 
section 6 the implications of this when comparing the results of 
the simulations with the observations. 

Fig. [6] shows the initial distribution of the density and Alfven 
speed va (isothermal sound speed) with some field lines of the 
magnetic field superimposed. Table [T] gives the Alfven speed at 
four locations for the two initial conditions considered (see sec- 
tion 5). Note that va increases as a function of height (away from 
the fluxrope), but varies less inside the loop. 

The flux rope is made to erupt by invoking an artificial force 
that acts on the flux ro pe plasma elements du ring the simulation, 
similar to the study of Pomoell et al. ( |2008| ). The Cartesian grid 
of the simulation consists of 300 x 300 cells in the x x j = [-0.2 
Ro, 0.2 Rq] X [1.0 Rq, 1.5 Rq] simulation domain. The boundary 
values are kept fixed during the simulation. 
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Fig. 6. The initial state of the simulation: density (left half) and 
Alfven speed (right half). The white lines depict magnetic field 
lines. The units in the color bar are for the density (left values 
in the color bar) 1.67 x 10"^^ gcm"^ and for the Alfven speed 
(right values in the color bar) 143 kms"^ the isothermal sound 
speed. Note the clipping of the color bar; a brown color indicates 
values larger than the color bar maximum. 



(a) 



(b) 



(c) 




Fig. 7. The density (left half) and speed (right half) at three dif- 
ferent times, t = 220 (a), t = 270 (b) and t = 330 (c) seconds. 
The white lines depict the magnetic field lines. The units in the 
color bar are for the density 1 .67x 10"^^ g cm"^ and for the speed 
143 km s"^ Note that both quantities are plotted using the same 
color bar. 



5. Results of the simulations 

Fig. [7] shows the mass density and speed at three diff'erent times, 
t = 220, 270 and 330 seconds. The dynamics of the eruption is as 
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Fig. 8. Left: Plot of the height- time data of the flux rope (solid line) and shock front (broken line) for the three simulation runs. Red 
corresponds to run A, green to run B and blue to run C. Right: Plot of the corresponding flux rope and shock speed. 



follows: As the flux rope starts to rise, a perturbation is formed 
around the flux rope. Due to the gradient in the Alfven speed 
and the increasing speed of the flux rope, the wave steepens 
to a shock ahead of the flux rope. However, the strength of the 
shock remains weak in the area below the loop. When the shock 
reaches the dense loop, it strengthens and slows down quickly 
due to the low Alfven speed in the loop (panel a). The erupting 
filament continues to push the loop structure ahead of it, acting 
as the driver of the shock. Thus, the speed of the shock is roughly 
that of the displaced loop structure when propagating in the re- 
gion of low va- As the filament decelerates, the displaced loop 
and shock escape from the filament (panel b). When reaching 
the region of higher y^, the speed of the shock increases with the 
increasing y^, and escapes from the propagating loop structure 
(panel c). 

To quantify the description above, we plot in Fig. [8] the 
height-time profile of the filament and the shock. We consider 
three diff'erent runs with slightly diff'erent parameters. Run B 
(blue color coding) and C (green color coding) are identical but 
for the magnitude of the artificial accelerating force, which is 
larger for run C. Run A has the same magnitude of the accelerat- 
ing force as run C, but the density contrast between the loop and 
ambient corona as well as the strength of the magnetic field is 
lower for run A. The speeds of the flux rope and shock are also 
shown in Fig. [8] 

In all three cases, the shock is strong while it propagates in 
the loop, with the compression ratio remaining between 3-3.5. 
However, the behaviour after the shock exits the loop depends on 
the simulation run. For run A, the shock only slightly decreases 
in strength, while for runs B and C the compression ratio quickly 
decreases to approximately 2.2. Interestingly, the shock strength 
remains at this value until the end of the simulation. 

In Fig. [9] we plot the radio emission produced by the shock 
assuming that the emission is produced immediately in front of 
the leading edge shock. The emission lane shows a frequency 
drop in accordance with the exponentially decreasing density. 




50 100 150 200 250 300 350 400 

Time, seconds 

Fig. 9. Radio track of the three simulation runs. The size of the 
marker indicates the compression ratio of the shock. 



Additionally, we indicate the compression ratio of the shock by 
the size of the marker. 

In addition to performing simulations with one dense loop, 
we also considered simulations with two loops. Adding another 
dense loop above the existing loop does not alter the overall dy- 
namics of the eruption. However, the shock structure becomes 
very complex as the shock enters the second loop. We refrain 
from showing the results, but note that a similar plot as Fig. |9] 
in the two-loop case shows two distinct tracks, corresponding to 
the shock propagating in the dense loops. In between the loops 
the frequency drops significantly. 
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6. Discussion and conclusions 

Observations on 13 May 2001 show that: 

- The speed of the white-light, plane-of-the-sky CME 
(430 km s"^) was similar to the projected speed of the erupt- 
ing filament (380 km s'^) and the EUV "blob" (450 km s'^). 
The CME most probably consisted of the ejected filament, 
and the filament propagated in the wake of the soft X-ray 
loop-like front 

- The projected speed of the soft X-ray loop-like front 
(650kms~^) and the estimated speed of the type II burst 
source (> 850 km s~^) were higher than the CME and EUV 
structure speeds 

- The estimated soft X-ray loop densities agree with the 
plasma densities of the fragmented metric type II burst. 
Also magnetic field strengths are in agreement with a super- 
Alfvenic shock at speed « 850 km s~^ near the soft X-ray 
loop heights 

- The fragmented type II burst changed appearance to a "reg- 
ular" type II burst at a time when the propagating soft X-ray 
loop-like front approached the outer boundary of the active 
region 

- The type II burst disappeared abruptly, without any weaken- 
ing of the emission. 

The MHD simulations that we have performed for a flux rope 
accelerated in a corona with dense loops show that: 

- The CME eruption drives a shock wave ahead of it through 
the corona 

- The shock wave propagates faster through the corona than 
its driver 

- In the dense loops, the shock becomes significantly stronger 

- Under right conditions, the shock strength decreases rapidly 
as the shock exits the loop. 

It is consistent with the simulation model that the high speed 
of the loop-like soft X-ray front actually corresponds to a strong 
shock wave propagating through the loop. When outside the ac- 
tive region, the changing geometry no longer favours shock ac- 
celeration and the type II burst dies out. Together with the as- 
sumption that the shock wave generates conditions that produce 
type II bursts, our MHD model produces burst lanes that cor- 
respond qualitatively to those observed in the event. However, 
in order for the radio signature to become fragmented as is ob- 
served, the conditions for plasma emission have to be some- 
how more favourable inside the loop than in the interloop area. 
The obvious assumption, consistent with our simple simulation 
model, is that the shock strength decreases significantly in the 
space between the denser loops. This requires specific condi- 
tions, such as a very large difference in the Alfven speed between 
the loop and interloop area, or a rapidly increasing Alfven speed 
in the interloop area. 

Our simplified two-dimensional model of coronal loops gen- 
erates the favourable conditions for shock formation as it neces- 
sarily yields weak magnetic fields in the dense loops to achieve 
lateral pressure balance. This is obviously too restrictive an as- 
sumption, as real coronal loops involve both increased densi- 
ties and magnetic fields. Thus, the strength of the shock in our 
model does not necessarily carry over to a more realistic model 
of the corona. However, the overlying coronal structure may ac- 
tually be quite fragmented containing loops with different densi- 
ties and magnetic fields, and our simulation model indicates that 
those with the lowest Alfven speed produce the strongest shock. 
Note also that the strength of the shock is not the only factor 



affecting the radio emission from the shock. Other factors, such 
as the density of the electrons in the loop, are important for the 
plasma emission mechanism as well. Thus, the densest loops in 
the corona overlying the active region seem most probable sites 
of plasma emission. 

It is also possible and even plausible that the three- 
dimensional structure of the corona plays a significant role. For 
instance, the different fragments of the bursts could originate 
from different loops that are not on top of each other. However, 
even in this case the interloop area must have specific conditions 
but not necessarily related to the strength of the shock, in order 
to cause the burst to suddenly stop emitting. 

In conclusion, using results from multi-wavelength data 
analysis and numerical MHD simulations, we have established 
that the unusual metric type II radio burst on 13 May 2001 can be 
explained by a model, where a coronal shock driven by a mass 
ejection passes through a system of dense loops overlying the 
active region, where the ejection originates. 
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